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1.  Introduction 


1.1  Shape  Memory  Alloys  (SMAs) 

SMAs  are  a  class  of  materials  that  undergo  a  martensitic  transformation  under  a 
mechanical  or  thermal  stress.  The  classic  example  is  equiatomic  nickel-titanium, 
commonly  called  nitinol  or  NiTi,  which  was  developed  at  the  Naval  Ordinance  Lab 
in  1959.'  At  temperatures  exceeding  100  °C  NiTi  exists  in  the  cubic  austenitic 
structure;  when  cooled  to  room  temperature  the  material  undergoes  a  diffusionless 
transformation  to  a  body-centered  tetragonal  structure.  The  material  can  be  heated 
to  the  austenitic  state  and  cooled  to  the  martensitic  state,  rendering  a  shape  change 
as  it  transforms  between  phases  (Fig.  1). 


Austenite 


Martensite  Martensite 

Fig.  1  Austenite  to  martensite  transformation 

In  addition  to  the  shape  memory  effect,  this  phase  transformation  allows  for  unique 
material  properties  that  make  SMAs  attractive  options  for  specific  applications. 
Two  such  properties  are  superelasticity  and  the  ability  to  perform  mechanical  work. 
Superelasticity  or  “pseudo-elasticity”  refers  to  the  ability  of  the  SMA  to  recover 
large  amounts  of  strain.  Generally,  NiTi  can  recover  up  to  6%  strain  without  plastic 
deformation.^  This  allows  the  material  to  have  quite  a  high  work  output.  The 
actuation  action  that  occurs  between  the  austenite  and  martensite  can  be  used  to 
store  mechanical  energy.  Together,  these  properties  give  SMAs  a  distinct  advantage 
over  other  materials  in  specific  applications. 

There  are  several  key  temperatures  that  help  define  the  phase  transformation: 
martensitic  start  (Ms),  martensitic  finish  (Mf),  austenitic  start  (As),  and  austenitic 
finish  (Af).  The  thermal  hysteresis  shown  in  Fig.  2  illustrates  the  general  regimes 
where  these  temperatures  are  defined.  If  the  material  starts  in  the  austenitic  phase 
and  is  cooled,  it  will  begin  to  form  when  the  temperature  reaches  Ms  and  continue 
until  the  Mf  temperature  when  the  material  is  completely  transformed.  At  this  point 
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it  can  be  heated  and  at  the  As  temperature  the  cubie  strueture  will  nueleate  and  form 
until  Af  is  reaehed. 


Fig.  2  Thermal  hysteresis  of  an  SMA 

This  thermal  hysteresis  defines  what  regime  the  SMA  ean  funetion  within.  If  the 
desired  applieation  requires  a  different  temperature  regime,  the  SMA  needs  to  be 
tuned  to  match.  One  method  of  tuning  SMAs  involves  ereating  nanopreeipitates  in 
the  matrix.  Nanopreeipiatates  provide  several  advantages,  they  suppress  rateheting 
during  thermal  eyeling,  redueing  the  possibility  for  failure  due  to  ereep. 
Additionally,  they  shift  transformation  temperatures  to  a  higher  regime  and  do  not 
interfere  with  the  transformation  proeess.^ 

The  transformation  temperatures  are  eompositionally  sensitive,  as  demonstrated  in 
Fig.  3.  Here,  the  dependenee  of  Ms  is  plotted  as  a  function  of  nickel  (Ni)  eontent 
and  it  ean  be  seen  that  as  Ni  shifts  from  50%  to  51%,  Ms  is  redueed  from 
approximately  75  °C  to  -150  °C.  It  would  seem  logical  to  use  titanium  (Ti)-rieh 
eompositions  instead;  however,  the  Ni  provides  stability,  strength,  and  so  forth. 
One  way  to  eombat  this  issue  is  to  introduee  preeipitates  by  aging  the  material.  NiTi 
can  be  aged  to  ereate  NfiTis  preeipitates,  whieh  generally  elevate  transformation 
temperatures."^  However,  to  push  SMAs  to  higher  temperatures  a  ternary  system  ean 
be  ereated.  The  addition  of  the  third  element  results  in  Ni-rieh  preeipitates  that  leave 
an  Ni-depleted  matrix.  A  matrix  with  an  Ni  defieit  will  have  higher  transformation 
temperatures. 
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Fig.  3  Transformation  temperatnre  vs.  Ni/Ti  ratio 

Due  to  the  nature  of  the  shape  memory  effect  and  superelacticity,  SMAs  are  a 
unique  class  of  materials  that  offer  a  unique  solution.  In  the  medical  field,  SMAs 
are  particularly  attractive  due  to  their  biocompatibility.  Perhaps  the  best  well- 
known  application  is  the  use  of  NiTi  wires  for  medical  stents.^  The  wire  is  tuned  to 
the  temperatures  of  the  human  body — it  can  be  inserted  in  a  smaller  shape,  and  then 
expanded  once  it  has  reached  the  desired  location.  The  constant  temperature  of  the 
body  ensures  that  the  stent  maintains  it  shape,  due  to  a  constant  stress  on  the 
material.  They  have  also  been  used  as  bone  staples  and  porous  bone  implants.^ 
Other  common  applications  include  orthodonture  wires  and  eyeglass  frames,  where 
a  constant  stress  is  put  to  work.  Although  SMAs  are  a  big  player  in  the  medical 
field,  they  also  find  applications  in  other  fields  such  as  aerospace^  and  energy.^ 


1.2  Jet  Turbine  Engines 

Increasing  the  efficiency  of  jet  turbine  engines  has  the  potential  for  great  economic 
and  environmental  improvements.  The  typical  jet  engine  consists  of  5  main 
components:  fan,  compressor,  combustor,  turbine,  and  nozzle.  Air  is  brought  into 
the  system  through  the  fan,  which  directs  some  of  the  flow  to  the  core  and  the  rest 
outside  to  act  as  a  coolant  and  noise  suppressant.  Air  directed  to  the  core  moves  to 
the  compressor,  where  its  volume  is  reduced  to  increase  pressure  and  potential 
energy.  The  condensed  air  is  then  fed  into  the  hot  section  of  the  engine  starting  with 
the  combustion  chamber.  Here,  the  high-energy  air  is  mixed  with  fuel  and  ignited; 
the  hotter  air  will  expand  and  is  therefore  storing  more  energy.  As  the  hot  air  moves 
through  the  next  section,  the  turbine  blades  rotate  and  move  blades  connected  via  a 
rod  in  the  fan  and  compressor.  Finally,  the  hot  air  exits  the  engine  and  mixes  with 
the  cold  air  from  the  initial  intake,  providing  thrust.^ 
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There  are  many  aspects  that  can  be  modified  to  improve  the  engine’s  efficiency. 
The  nozzles  injecting  fuel  can  have  different  shapes  or  hole  geometry  or  the  blade 
clearance  throughout  the  engine  can  be  modified,  for  example.  Another  promising 
approach  involves  articulating  the  blades  in  the  turbine  wheel.  Currently,  there  is 
no  way  to  articulate  the  blades  so  the  angle  of  pitch  is  one  that  is  broadly  beneficial 
for  ascension,  cursing,  and  dissension.  If  the  blades  could  be  articulated  such  that 
for  each  of  the  3  stages  they  are  in  a  different — but  optimal — ^position,  the  fuel  will 
be  used  more  efficiently. 

To  introduce  an  actuation  system  into  the  turbine  there  are  several  material 
properties  that  need  to  be  considered.  The  environment  will  be  high  temperature 
and  pressure,  and  will  require  a  high  power  density,  long  lifetime,  low  cost,  and 
lightweight  material.  Traditional  actuation  systems  such  as  hydraulic  and 
pneumatic  could  provide  the  desired  power  output,  but  are  generally  too  heavy  and 
bulky  to  incorporate.  Other  material  options  include  piezoelectric,  electromotive 
force,  or  smart  materials.  As  a  subset  of  smart  materials,  SMAs  are  a  lightweight, 
low-cost,  and  power-dense  option.  The  challenge  lies  in  bringing  them  to  function 
in  high  temperatures  and  pressure  with  a  reliable  lifetime. 

2.  Methods,  Assumptions,  and  Procedures 

This  work  was  performed  using  phase  field  analysis  on  the  Ohio  Oakley 
Supercomputer  and  Army  Excalibur  Supercomputer.  Phase  field  codes  are  explicit, 
numerical  solutions  that  replace  interface  boundary  conditions  with  a  partial 
differential  equation.  This  phase  field  code  uses  a  forth  order  Landau  polynomial 
(Eq.  1)  for  the  asymmetric  transformation  pathway  between  the  austenite  and 
martensite  phases  of  the  material.  The  free  energy  equation  is  calculated  with 
elasticity,  chemical  and  interfacial  information  as  detailed  in  the  (Eq.  2).  This 
simulation  reads  information  in  a  VUMAT  file,  which  provides  Abaqus  with  user- 
specified  material  parameters,  here  specific  to  NiTi.  Additionally  the  simulation  is 
provided  with  the  slip  systems  allowing  for  the  possibility  of  plastic  deformation. 
Einally,  the  user  can  define  the  mesh  used  including  desired  size,  shape,  and 
orientation  of  a  precipitate. 
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Eq.  1 ;  Landau  polynomial 


fLp{yt>  d}  —  ^ 


+  q 


+ 


(1) 


Eq.  2:  Eree  energy  equation 

Nt  Nt 

Wi,=-E^CE+f,p(,Vt,e)+-Y^Vvt 

t=l  U=1 


y  KtuVvt  .  (2) 


The  simulation  in  this  work  took  advantage  of  the  fast  eomputational  time  of  2-D 
geometries.  This  representation  ean  be  derived  from  a  3-D  model.  In  3-D  a  eube 
represents  the  austenitie  state  and  a  brick  represents  the  martensitic  state.  The  brick 
comes  from  stretching  the  cube  in  x,  y,  or  z  (Eig.  4).  Similarly,  in  2-D  austenite  can 
be  represented  by  a  square  and  martensite  can  be  represented  by  a  rectangle, 
stretched  in  x  or  y.  Boundary  conditions  include  a  fixed  origin  and  rolling  condition 
on  the  bottom  left  node.  The  material  starts  at  350  K  and  is  immediately  step- 
quenched  to  150  K  at  the  first  time  step,  then  held  for  the  remaining  time.  This 
allows  the  transformation  between  austenite  and  martensite  to  be  as 
computationally  efficient  as  possible.  Nucleation  sites  determine  where  the 
martensite  will  begin  to  form. 


Fig.  4  3-D  representation  of  anstenite  to  martensite 
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In  terms  of  precipitate,  orientation,  size  and  shape  were  values  of  interest  to  study. 
The  goal  is  to  look  at  the  effect  that  these  3  parameters  have  on  microstructure, 
transformation  temperatures,  and  actuation  properties. 

The  first  study,  detailed  in  this  report,  is  a  qualitative  study  of  the  effect  of 
orientation  on  microstructure.  An  elliptical  precipitate  with  a  long  axis,  a  =  0.25, 
short  axis  b  =  0.10  (according  to  Eq.  3),  was  rotated  from  0°  to  45°  in  15° 
increments.  Images  of  the  final  microstructure  were  visualized  in  Abaqus  CAE, 
version  6.14.  A  case  without  a  precipitate  was  also  run.  Each  case  was  run  once. 


Eq.  3;  The  form  of  an  ellipse,  long  axis  a,  short  axis  b, 


— +  —  =  1 


(3) 


Because  of  the  output  Abaqus  can  handle,  a  single  variable  can  be  plotted  at  a  time. 
This  means  that  if  “Variant  1”  is  plotted  in  red,  the  blue  is  simply  the  absence  of 
Variant  1  (as  opposed  to  directly  being  Variant  2).  Likewise,  if  “Variant  2”  is 
plotted,  elements  identifying  as  Variant  2  will  be  red  and  elements  not  identifying 
as  Variant  2  will  be  blue.  For  the  purposes  of  this  qualitative  study,  it  can  be 
assumed  that  red  is  Variant  1  and  blue  is  Variant  2  for  simplicity. 


In  the  final  step,  the  general  structure  of  the  martensite  can  be  thought  of  as 
“compatible”  or  “incompatible”.  Here  the  difference  between  compatible  and 
incompatible  is  the  alignment  of  the  precipitate  and  martensite  plates.  If  there  is  a 
good  alignment,  stress  will  be  reduced  locally  creating  a  compatible  situation.  If 
there  is  not  a  good  alignment — ^which  could  lead  to  plates  stretching  or  shrinking 
to  accommodate — the  precipitate  will  experience  a  local  stress,  in  an  incompatible 
situation.  The  benefits  or  detriments  of  stress  will  be  addressed  later,  but  are  helpful 
for  distinguishing  the  difference  between  compatible  and  incompatible. 


3.  Results  and  Discussion 


The  initial  study  investigated  the  distribution  of  martensite  around  the  precipitate 
as  a  function  of  orientation.  As  previously  detailed,  the  precipitate  has  a  long  axis 
and  a  short  axis.  The  patterns  formed  relative  to  each  axis  are  discussed  with  regard 
to  the  “interruption”  of  one  variant  by  the  other. 

In  the  case  with  no  precipitate,  the  final  step  shows  the  plot  of  Variant  1  in  red.  The 
blue  areas  are  assumed  to  be  Variant  2.  In  this  case  the  2  variants  form  an  alternating 
striated  pattern,  which  is  expected  as  the  material  experiences  twinning  in  the 
martensite.  Generally,  the  pattern  is  continuous  (Fig.  5). 
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Fig.  5  Phase  field  simulation  without  a  precipitate 

Now  we  add  in  the  precipitate  at  0°,  the  material  is  cooled  and  the  striated  pattern 
form  again.  This  time,  there  appears  to  be  some  disruption  in  the  line  orthogonal  to 
the  long  axis  of  the  precipitate.  Here,  variants  span  approximately  half  the  length 
expected,  only  to  be  interrupted  by  the  alternate  variant  type.  Although  the 
interrupted  striations  are  most  notable  in  this  direction,  there  is  interruption  in  the 
direction  orthogonal  to  the  short  axis  as  well  (Fig.  6). 


Fig.  6  Phase  field  simulation  with  precipitate  oriented  at  0°  to  the  horizontal  axis 

Next,  the  precipitate  is  rotated  to  15°.  Here,  the  interrupted  variants  have  a  much 
smaller  width,  particularly  closer  to  the  precipitate.  It  seems  as  though  the 
interruption  does  not  extend  as  far  through  the  material  as  in  the  0°  case,  but 
remains  more  local  (Fig.  7). 
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Fig.  7  Phase  field  simulation  with  precipitate  oriented  at  15°  to  the  horizontal  axis 

As  the  precipitate  is  moved  to  the  30°  position,  the  interruption  in  martensite 
variants  in  the  direction  parallel  to  the  long  axis  has  almost  completely  disappeared. 
Small  interruptions  exist  very  locally  around  the  ends  of  the  precipitate,  where  the 
plates  are  again  smaller  than  the  0°  case  (Fig.  8). 


Fig.  8  Phase  field  simulation  with  precipitate  oriented  at  30°  to  the  horizontal  axis 

Finally,  the  precipitate  is  rotated  to  45°  where  it  appears  to  be  aligned  with  the 
martensite  plates  that  form.  Although  there  appears  to  be  little  interruption  on  the 
long  axis,  there  is  minor  interruption  parallel  to  the  short  axis  on  a  very  local  scale. 
Overall,  the  interruption  does  not  extend  past  2  plates,  keeping  it  relatively  local 
(Fig.  9). 
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Fig.  9  Phase  field  simulation  with  precipitate  oriented  at  45°  to  the  horizontal  axis 

The  interruption  is  of  interest  because  the  existence  of  an  interface  makes  stress 
concentration  development  or  residual  strain  more  likely.  If  the  edge  of  the 
precipitate  is  in  an  incompatible  alignment  with  the  variant  pattern  in  the  matrix, 
there  is  likely  stress  at  the  interface  between  precipitate  and  matrix.  Although 
plastic  deformation  in  the  form  of  residual  strain  is  not  desirable,  the  ability  to 
control  what  type  of  variant  forms  where  is  desirable.  If  the  local  stress 
concentrations  around  the  precipitate  can  dictate  which  type  of  variant  forms,  it  can 
select  for  the  more  compatible  variant  during  the  next  thermal  cycle.  This 
essentially  allows  the  precipitate  to  become  a  nucleation  site.  A  controlled 
microstructure  means  a  more  stable  material  and  in  turn  a  reliable,  long-lasting 
option  for  high- temperature  actuation  applications. 

4.  Conclusions 


Creating  an  actuation  system  that  can  move  turbine  blades  during  flight  of  an 
aircraft  requires  very  specific  material  properties.  SMAs  are  an  excellent  option  as 
they  have  a  high  energy  density;  however,  they  are  currently  lacking  the  high- 
temperature  stability  necessary  for  operation  in  extreme  thermal  conditions.  To 
push  SMAs  to  higher  temperatures  nanoprecipitates  are  formed,  allowing  the 
martensitic  phase  transformation  to  stabilize  over  many  thermal  cycles  while 
increasing  the  transformation  temperature  range.  How  these  precipitates  work  is 
not  well  understood,  but  harnessing  their  features  aids  the  design  of  a  new  high- 
temperature  SMA.  This  work  used  phase  field  analysis  to  show  that  the  orientation 
of  the  precipitate  is  of  special  interest  due  to  the  effect  on  the  microstructure.  A 
precipitate  oriented  with  the  martensitic  variants  shows  little  potential  for  residual 
strain,  which  may  lead  to  failure.  If  the  precipitates  are  at  a  desired  orientation,  the 
phase  transformation  may  be  able  to  occur  safely  at  much  higher  temperatures  than 
previously  thought. 
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5.  Recommendations  and  Future  Work 


Before  SMAs  ean  be  implemented  into  high-temperature  actuation  systems,  they 
need  cyclic  stability.  Through  the  use  of  this  phase  field  analysis,  the  local  stress 
around  a  precipitate  can  be  used  to  create  a  stable  material.  Quantitatively,  there  is 
still  much  work  to  do  in  determining  local  stresses.  Once  the  optimal  precipitate’s 
size,  shape,  and  orientation  have  been  selected  due  to  their  effect  on  the 
microstructure,  transformation  temperatures,  and  actuation  properties  through  this 
phase  field  analysis;  then  the  information  can  be  fed  into  a  larger  scale  model.  Due 
to  the  application  of  creating  an  SMA  for  turbine  blade  actuation,  predictions 
beyond  the  scale  of  phase  field  calculations  is  necessary.  Microstructural  finite- 
element  simulations  perform  calculations  on  a  larger  scale.  This  model  can  be 
informed  by  the  phase  field  output  as  well  as  data  collected  from  experimental  work 
done  at  National  Aeronautics  and  Space  Administration’s  Glenn  Research  Center 
and  The  Ohio  State  University.  Together  this  information  will  provide  a  material 
that  can  perform  work  at  high  temperatures,  which  can  be  fed  into  the 
microstructural  simulation  that  will  speak  to  the  feasibility  of  implementing  SMAs 
as  actuators  in  jet  turbine  engines. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


2-D 

2-dimensional 

3-D 

3 -dimensional 

Af 

austenitic  finish 

As 

austenitic  start 

APG 

Aberdeen  Proving  Ground 

HPC 

high  performance  computing 

Mf 

martensitic  finish 

Ms 

martensitic  start 

Ni 

nickel 

osu 

The  Ohio  State  University 

SMA 

shape  memory  alloy 

Ti 

titanium 
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